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Abstract 

The Cauchy problem for the Gross-Pitaevsky equation with quadratic nonlocal nonlinearity is 
reduced to a similar problem for the correspondent linear equation. The relation between symmetry 
operators of the linear and nonlinear Gross-Pitaevsky equations is considered. 



Introduction 

Recent advances in formation of the Bose-Einstein condensates (BECs) of alkaline metal vapors 
[1-3] have stimulated study of the theoretical models describing behavior of nonlinear systems in 
external fields. In the BEC models the local and nonlocal Gross-Pitaevskii equations (GPEs) 
(which are also called the Hartree-type equations in the mathematical literature) are widely used. 
Besides the BEC theory, the nonlocal GPE serves as a basic equation in the models of quantum 
many-particle systems, nonlinear optics, collective excitations in molecular chains, etc. 

Let us write down the nonlocal Gross-Pitaevskii equation as 

{~ihdt + HUt)}'9{x,t) = {-ihdt + n{t) + >cv{t,'f{t))}'f{x,t) ^ 0, (i) 

^{x,t)eL2m), V{t,^it))^ f dy^*{y,t)Viz,w,t)^iy,t). (2) 



Here the linear operators 7i (i) = Ti-iz, t) and V{z, w, t) are the Weyl-ordered functions [4] of time 
t and of noncommuting operators 

z = {p,x) — {—ihd/dx,x), w — {—ihd/dy,y), x,y ^ M", 

with commutators 



[zk.Zj]^ = [wk,Wj]^ ^ ihJkj, [zfe,ii)j]_ = 0, fc,j = l,2n, (3) 
J = ll^fcj ||2nx2n IS & identity symplectic matrix J = \ ^ ) i I = Inxn is an identity 

V ^ /2rix2n 

n X n-matrix. 

In multidimensional space GPE ([T|) with variable coefficients of general form is non-integrable 
by known methods like, e.g., the Inverse Scattering Transform [5]. Therefore, analytical solutions 
of this equation can be constructed only approximately. An effective approach to construct such 
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solutions is provided by the method of semiclassical asymptotics. Thus for nonhnear self-consistent 
field operators, the theory of canonical operator with real phase was constructed for solution of 
the Cauchy problem in [6,7]; for spectral problems including singular potentials in [8,9] (see also 
[10-13]). Soliton-like solutions of the Hartree-type equation with some potentials of special form 
were constructed in [14]. A specific and attractive feature of the nonlocal GPE with nonlinearity 
presented in the equation only as a term under integral sign is that this equation can be referred to 
a class of nonlinear equations of mathematical physics which are close to linear ones in a sense [4] . 
Namely, among solutions there is a subset of solutions regularly depending on the nonlinearity 
parameter. Therefore, on the class of functions called trajectory-concentrated functions [15], 
the problem of semiclassical asymptotics construction for the nonlinear equation is reduced to 
an auxiliary problem of construction of asymptotic solutions for associated linear Schrodinger 
equations. Finally, the WKB-Maslov complex germ method [16, 17] was generalized for Eq. ([T]). 
In particular, formal solutions of the Cauchy problem asymptotic in formal small parameter h 
{h — > 0) was constructed accurate to 0{h^^^) where N is any natural number. The leading term 
of asymptotic solution of the spectral problem was found in [18]. 

Let us note that the semiclassical method, being approximate in essence, allows one to find 
exact solutions in some special cases. In [19, 20] the evolution operator for Eq. ([T]) with the 
quadratic nonlocal potential was found. In this case the linear operators t) and V{z, w, t) are 
quadratic in i, w: 

n{z,t) = ^{z,HUt)^) + {n^{t),z), (4) 

V{z,w,t) = ^{z,W,,{t)z) + {z,W,^,{t)w) + ^{w,W^Ut)w). (5) 

Here, Hzz{t), Wzz{t), Wzw{t), Www{t) are 2n x 2n-matrices, Hz{t) is a 2n-vector; (.,.) is an 
Euclidean scalar product of vectors: (p, af) = ^ Pj^jj Pi 2? € K", {z, w) = "j^ z,w £ R^". 

The Hartree-type equations ([1]), (j4|), ^ are integrated explicitly [19,20] and possess fairly 
rich symmetries, the study of which gives varied information about solutions of the equation. It 
also shows how to reduce multidimensional equations to one-dimensional equations, to construct 
classes of exact and approximate solutions, to investigate asymptotic of special classes of solutions, 
etc. Moreover, as long as Eqs. ([1]), ([4]), ([5]) have nonlocal nonlinearity, its symmetry is of a 
special interest in the symmetry analysis. The matter is that standard methods of symmetry 
analysis [21-25], developed basically for partial differential equations (PDEs), face the problems 
when structure of equation under the study differs from the PDE because, for example, there are 
no regular rules to choose an appropriate structure of symmetries for a non-differential equation. 
Equation ([l]) enable one to avoid this problem as symmetry of this equation is closely connected 
with the symmetry of linear equation corresponding to the nonlinear one. 

Here, we find symmetry operators of Eqs. ((Il), (|4]), ([SJin explicit form which, by definition, leave 
invariant the solution set of the equation and allow to generate new solutions from the known ones, 
(see, e.g., [26,27]). 

1 Operators in terms of an operator Caushy problem 

Consider the Cauchy problem for Eq. ([T]) 

■^{x,t,h)\t=o^-i{x), 7e§, 117(2^)11' = 1 (6) 

where § is a Schwartz space. 

The second-order Hamilton-Ehrenfest related to the Cauchy problem ^ ^ has the form [20] 

= Ji^zit) + [n^t) + K{Wzz{t) + Wz^iMz^}, 

(7) 

A*2 - J[nzz{t) + xWzz{t)]^^2 - Avt2[W.^(t) + KWz^{t)]J. 

Let us change the function 5'(af, t) by ^(a;, t) in Eq. ([1]) as 



^ {x, t) = exp { [S{t) + {P{t) ,x-X{t))]]'^{x-X{t),t), 
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(8) 



where S{t), P{t), X{t) are some differential functions to be determined. 
For ^{x, t) we have the following equation: 

{ - ihdt + Sit) + {P{t),x-x{t)) - {P{t),x{t)) + i(z$,H,.(t)z$> + {nAt),z^) + 

where we use the notations z$ = {—ihd/dx + P{t),x); w$ = {—ihd/dy + P{t),y). 

By changing variables x = u + X{t) and denoting <i>(u + X{t),t) = t), we obtain the 
following equation: 

- ihdt + ih{Jc{t), ^) + s{t) + {P{t),u) - {P{t)j{t)) + hz^u.num^u) + 

ou 2 

+ i(w$„,iy^^(i)u)$„)^l>(y,t)|$(u,t) =0, (10) 

where = z„ + Z{t), — Zy + Z{t), z„ = {—ihd/du, u), Zy = {—ihd/dy, y). 
The functions '^{u,t) are centered by construction, i.e., satisfy the condition 

^*{u,t)zu^u,t)du^Q. (11) 



Initial condition (jB]) for Eq. has the form 

|.(u,i)L^o =exp{ - '-[S{0) + {P{0),u)]yi{u + X{{)))- (12) 
Let the vector z = Z{t) = {P{t), X{t)) satisfies the equation 

z{t) - j{n,{t) + [n,,it) + HWz.it) + w,^it))]z{t)} (13) 

with the initial condition Z{0) = ('y{x)\z\"f{x)), and the function S{t) is determined by the relation 

t 

Sit)^ j[{P{t),k{t))-m]dt, (14) 



where 

?){t) = ^{Z{t), [nUi) + HW,,{t) + 2W,^{t) + W^^{t))]Z{t)) + 

+{n,{t), z{t)) + ii^Sp(w^^(t)A2). 

Here the matrix A2 of the order 2n x 2n satisfies the equation 

A2 = mut) + i(W,,{t)]^2 - A2[n,,{t) + ^wut)]J (15) 

and the initial condition 



A2(0) = -||(7(f)|{A%Azfc + AzfeA%}|7(f))||. (16) 



Then the function ^{u,t) is a solution of the linear associated equation 

{ - ihdt + {n,,{t) + KW,,)zu)]Hu,t) = (17) 



with the initial condition fT^ . 

Let A{u, t) — A{zu, t) be the operator whose Weyl symbol A{z, t) satisfies the relation 

-ihdt + ^{zu,{'Hzz{t) + >cW^^)zu),A{u,t)'\ =0 (18) 
and the initial condition 

A{u,t)\t^o = aiu), (19) 

where a{u) : § ^ § is an arbitrary operator. 

Then, thereby (fTB]) the operator A{u, t) is a symmetry operator of Eq. PT|) that maps a solution 
^{u,t) of Eq. p?)) into another solution of this equation. Respectively, the function determined 
by the condition 

$a(u, t) = —A{u, t)^{u, t), (20) 
aA 

where aA = 0)|| is also a solution of Eq. ^T7\ . 

At t = we have 

¥a(m, 0) = —a{u)<j){u) = (j)A{u). (21) 
aA 

Here [[(/(^(w)!! — 1, and we immediately obtain ||$yi(M, t)|| — 1. However, in general case the 
function $a(w, <) can correspond no solution of the original nonlinear equation as it is centered 
not for all A{u, t): 

j C(u, 0)z„1a(m, 0)d?2 ^ 0. (22) 

To find solutions of the original nonlinear equation which would correspond to $(w, t), let us 
introduce the notations 

^= '^*A{u,0)zyWAiu,0)du = / (t)\{u)zu<f)A{u)du (23) 



and 

Kt) = / '^Aiu,t)z^'^A{u,t)du. (24) 



We can verify immediately that if X{t) is a solution of the Cauchy problem 

m - Jin^.it) + ^^WUmit), X{t) = (b'^^^^ , A(0) Ao - (^V^ , (25) 
then the function 

^Aiu,t) = exp{^[5A(i) + (Ap(t),u + A„(i))]}$A(M + An(i),t), (26) 

where 

t 

Sx{t) = J {{Xp{t),L{t))-Sjx{t)}dt (27) 



and 

S)x{t) - i(A(t), [H.At) + HW^t) + 2W,Ut) + W^^mX{t)) + ^KSp{W^Ut)^2) 
is a solution of Eq. (fT7)) and satisfies the condition 

$^(u, t)z„$A(u, i)c?w = 0. (28) 
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Let us correlate according to ^ the functions ^a{u, t) = <i>(x + X^it), t) and i) = $(5? + 
X{t),t) with the functions ^'^(a;, i) and ^'(a?', i), respectively. Then we have the relation 

exp{ - 'j^[SA{t) + {PA{t),X-XA{t))]]^A{S,t) = 
^—C^^{USx{t) + {\p{t),X-XA{t) + \u{t))]]A{x-XA{t) + \u{t):t)x 

UA '^n J 
X exp { - ^-[S{t) + {P{t),x - X{t))]^-^{x + X{t) - XA{t) + \u{t), t) (29) 

from which it follows 

A{x,t) = ex^Y>[\[SA{t) + {PA{t).S - XA{t))]] X 
X — exp|^[5A(t) + {\p{t),x^ XA{t) + \u{t))]]A{x- XA{t) + \u{t),t) X 

xexp| ~ ■^[5(0 + (P(t),a;-X(t))]}«'(f + -XA(0 + A„(0,i). (30) 

The operator 

*A(x,i) l„i*(f,i) 
determined by pop is a symmetry operator of the original nonlinear equation. 

2 Symmetry operators obtained via auxiliary equation 

Let a : 5' ^ S* be an operator and 7(x) £ 5 is a function. Let us set for Eq. ([1]) the following 
Cauchy problems: '^{x,t,h)\t=a — ^{x), ^{x) £ V^, \\l{x)\\'^ — 1, and 'ifA{x,t,h)\t=o = a^{x) = 
7A(f),7A(f)e7'0, ||7a(x)|P = 1.^ 

Let us change the functions ^(a;, t) and '5/1 (a?, t) in Eq. H]) by the functions $(«, t) and ^^(m, t), 
respectively, according to formulas ([5]), (fT^ . HH), and Then the function $(u, <) satisfies the 
equation 

I _ i^ia, + (7^^^(t) + >^t^,,)z„)|$ ^ (31) 

and the initial condition ^{u,t)\^_^^ — exp | — — [S'(O) + (P(0), u)] |7(m + ^(0)), where the vector 

Z(i) is a solution of Eq. (fT5|) with the initial condition Z(0) = (7(x)|z|7(x)). Similarly, $a(u, 
is a solution of the equation 

{ - ihdt + (7^..(i) + ><VF,,)z„)}$a(??,0 = (32) 

with the initial condition I>a(u, ^)|(^g = exp | - [5*^(0) + (Pa(0), u)] |7a(u + -?yi(0)). 
Vector ZA{t) satisfies Eq. (|13p with the initial condition Za{0) ~ {'yA{x)\z\jA{x)) , and 

t 

SA{t) = j[{PA{t),XA{t))-^A{t)]dt, (33) 


where 

^A{t) = \{ZA{t), [HUi) + ^{W,,{t) + 2M^,,„(t) + W^^{t))]ZA{t)) + 

+ (H.(t),ZA(t)) + ii^Sp(M^„„(t)A2A). (34) 
Matrix (t) is a solution of Eq. with the initial condition 

A2a(0) = i||(7A(f,i)|{AzjAzfc + AzfcA%}|7A(x,i))||. (35) 

5 



The relation 7^(5^) = a'y{x) determines operator a{u) such that 0) = a{u)^{u,0). 

Consider the Cauchy problem for an operator A{u,t) — A{zu,t): 

- ihdt + ^{zu, {n,,{t) + xW,,)zu),A{u,t)\ = (36) 

with the initial condition ^ 

'A{u,t)\t=o = a{u). (37) 

In accordance with A{u, t) is a symmetry operator of Eq. PT|) and it maps solution $(-u, t) of 
Eq. (pij) into its solution (m, t) . 

If we correlate the functions I>a(u, t) = <^{x -\- X^(i), i) and $(m, t) — (^[x + X{t),t) with the 
functions '^A{x,t) and \E'(a;, i), respectively, according to formula (H]), we obtain 



[~ ^[SA{t) + {pAit),X - XAm}'i'Aix,t) ^Aix - XA{t),t) X 

{ - ^[5(0 + {P{t),x - X{t))]}^{x + X{t) - XA{t), t) (38) 



exp 

X exp 



or 



*^(f , t) = exp { [5a (0 + {PA{t). x-XA{t))]]A{x- XA{t),t) X 

xe^'p[-'j\S{t) + {P{t),x-X{t))]^^{x + X{t)-XA{t),t). (39) 

An operator 

^^(^E,^) = A„i^'(f,i) 

determined by relation p9p is also a symmetry operator of the initial nonlinear equation. 

3 Conclusion 

It should be noted that a direct calculation of symmetry operators for a nonlinear equation is, 
as a rule, difhcult because of the complexity of the determining equations [28]. In construction of 
symmetry operators, we have used the fact that the original nonlinear equation can be connected 
with an associated linear equation, and for quadratic operators of the form Q and ([5]) these 
associated equations coincide. In general, this is not true and therefore, analytical solutions of the 
GPE can be constructed only approximately. 

Many approximate methods are based on appropriate anzats representing the general element of 
a class of functions in which approximate solution is constructed. To such methods it is possible to 
relate, e.g., the method of collective variables (see, e.g., [29] and references therein), the Lagrangian 
method [30, 31], etc. 

It appears that the most effective method in various multidimensional problems of mathematical 
physics is the method of semiclassical approximations. Peculiarity that differs this method it 
from the usual method of expanding in power series in a small asymptotic parameter is that 
the asymptotic small parameter is included into the solution both regularly and singularly. This 
fact allows us, in particular, to construct space-localized solutions having an important physical 
meaning. This possibility is of a special importance in nonlinear problems where stable localized 
excitations (patterns), such as, for example, solitons are the object of research. The advantage 
of the semiclassical asymptotic method is that it allows to estimate accuracy of the constructed 
solution with a given power of the asymptotic parameter. 

Symmetry analysis provides another approach to construct analytical solutions of the nonlinear 
equations [21-25]. Taking into account invariant properties of the equation under consideration, 
one can find classes of special solutions which can serve as a prototype of anzatses for classes of 
exact solutions. 

However, variable coefficients, as a rule, reduce the equation symmetry and even can exclude it 
completely that reduces the symmetry analysis in such cases. Comparison of these two approaches 
results in a problem of consideration of approximate symmetry for the GPE in the formalism of 
semiclassical approximation method (see, e.g., [32]). 
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